


AIAA Paper 
No. 73-701 


Tore) OF OLIN THE ARCTIC 


by 
DAVID) P, HOULT 
Massachusetts Institute of Technology 


Cambridge, Massachusetts 








First publication rights reserved by American Institute of Aeronautics and Astronautics. 
1290 Avenue of the Americas, New York, N. Y. 10019. Abstracts may be published without 
permission if credit is given to author and to AIAA. (Price: ALAA Member $1.50. Nonmember $2.00). 


AlIAA Library, 750 3rd Avenue, 


 POLARPAM 





New York, New York 


AIAA 6th Fluid and Plasma 
Dynamics Conference 


PALM SPRINGS, CALIFORNIA / JULY 16-18, 1973 


Note: This paper available at AIAA New York office for six months; 
thereafter, photoprint copies are available at photocopy prices from 


10017 


9°S99°GG°Z0G: weg 


NOH 


i ee 


un 


M4 3 


if eae 


os 





\ 





THE SPREAD OF OIL IN THE ARCTIC 


David P. Hoult 
Associate Professor 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 


Abstract 


The problem of oil spreading over Arctic ice 
is studied. The main topics considered are, first, 
what is the rate of spread of oil released froma 
point on the surface? Second, how large does the 
area of the oil slick become before it is caught 
in the pockets of the ice surface? The random 
character of the ice surface controls the final size 
of the oil slick and its rate of spread. The theory 
of the spreading speed essentially balances a drag 
force against a gravitational driving force. The 
spreading theory contains one emperically deter- 
mined constant. The experiments were performed in 
a small tank in the Fluid Mechanics Laboratory at 
Massachusetts Institute of Technology, using random 
surfaces constructed from plywood, gravel, and 
epoxy. The spread and final size of the spill, as 
observed, UES in agreement with the predictions of 
the theory - In addition, the results were 
compared with larger spills made in the high Arctic 
during two Pat oNey aoe carried out by the U.S. 
Coast Guard fe - It is found that,properly 
scaled, the field data are in good agreement with 
both the laboratory data and the theory. These 
results serve to provide an understanding of the 
potential for oil pollution when oil is spilled in 
the high Arctic on either tundra or pack ice with 
no leads. A qualitative argument is given to show 
when the results can be applied to a surface of 
modest slope. 


I. Introduction 


In 1968, oil was discovered on the North Slope 
of Alaska. Since the world supplies of oil are 
dwindling and the prices rising, this oil must be 
used if the present standard of living is to be 
maintained at a reasonable cost. This means that 
the crude oil must be shipped to a place where it 
can be refined, and whether supertanker or pipeline 
is used, the oil will have to travel over the rough 
surfaces of the Arctic. A small fraction of all 
oil shipped is spilled. In order to deal with this 
environmental problem it appears necessary to gain 
an understanding of how oil spreads on or under 
ice. This understanding must come in two parts: 
first, the process of spreading must be understood 
and the time for spreading determined; and second, 
the final size of the oil pool must be determined. 


In developing the spreading law, proper 
account must be taken of the characteristics of 
the rough surface. In order to solve the second 
problem, one must take into account the character- 
istics of the ice surface, namely its roughness and 
the random nature of this roughness. If this is 
done, the answer to the problem will be determined 
ina stapyetical way. This problem was studied by 
Patureau The purpose of this paper is to 
correlate results of experiments on a laboratory 
scale model and comparing these results with field 
data. 


The experimental apparatus and procedure will 
be discussed in Section II. Since the results of 
Patureau are expressed in terms of parameters 
derived from the power spectrum of the surface, 
this section also describes the apparatus and 
procedure used to determine the power spectrum of 
the experimental surface. In Section III, the 
details of Patureau's final size theory and the 
spreading theory derived from the experimental 
results will be discussed, along with the assump- 
tions and limitations. Section IV will present the 
results in the form of correlations between theory 
and experiment, an explanation of the derivation of 
the spreading theory, sources of error, and the 
conclusions of the study. 

II. Experiments 

In planning the experiments, the characteris- 
tics of Arctic ice were first studied. Earlier 
work had indicated that it was desirable to con- 
struct an apparatus thay could dump a known volume 
of ‘oil’ at a fast rate This turned out later 
to not be the case. 


To obtain dynamic and geometric similarity, 
certain scaling parameters must have suitable 


values. These are: 
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4 Sey, 1/80 above 
1/20 below ice 
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number of pockets >> 1 


A 


density ratio <.2 (under only) 


With the exception of Re, the conditions are all 
well met. In the case of the Reynolds number, a 
desire to use fluids that could be used easily, 
would not destroy plexiglass and would provide 
agreement with the other parameters led to a relaxa- 
tion of the restriction. Water for the over ice 
experiments provides an Re of a little more than 

400 and kerosene for the under ice case provides 

an Re of about 100. This is believed to provide 
sufficient turbulence to satisfy the theory. 


A schematic of the spreading experiment 
apparatus is shown in Fig. (1). This figure should 
be referred to if any questions arise. 


The rough surface itself was built by epoxying 
various sizes of rocks to a plywood board in a 
random fashion. The surface was then painted white 
with black concentric circles. When used in con- 
junction with a dyed fluid, the spread is easily 
observed and photographed. The size of the surface, 
and therefore the size of the tank,is 4' x 4‘. 

This size tank offers a large enough surface to be 
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statistically significant and yet small enough to 
be easily handled. 


In order to scale the laboratory model to a 
real spill, one should consider that the Arctic 
roughness height is on the order of 1' . The rough- 
ness height of the experimental surface is about 
1/10 inch. This means that the 4' x 4' surface 
corresponds to about 500' x 500' in the Arctic. 


Tha tank is made of plexiglass to allow the 
spread to be easily observed and photographed. The 
frame was designed to allow for rotation of the 
tank when conducting under ice experiments. This 
was necessary because if the tank were filled ina 
horizontal position the pockets in the surface would 
trap air and make the experimental results meaning- 
less. 


To record the progress of each experiment a 
Bolex 16mm camera with a 10mm lens and a variable 
time base was used. Since this lens needs a focal 
length of around nine feet to get a large enough 
field to cover the apparatus, it was decided to use 
a mirror. The camera was then run at a known speed, 
usually 20 fps, to record time. The developed film 
was then examined frame by frame to determine radius 
versus time and the final size. 


To measure the surface a frame was made (see 
Fig. 2) and the elevation of the surface was 
measured. The variance and the power spectrum were 
then calculated by a computer program. 


To test the isotropy assumption, traces were 
taken in perpendicular directions and then their 
power spectra compared. This is explained in 
Section IV. 


The frame affects the spectrum only in the 
wavelength range greater than ten inches. The 
imperfections in the frame has no effect on the 
results. 


For the over ice experiments, the fluid used 
was water dyed with food coloring. Soap was added 
to reduce the surface tension and the 'oil' was 
placed in the cannister. The apparatus was then 
leveled very carefully to prevent the fluid from 
running off to one side. Problems were encountered 
in this area. Once the experiment was ready to 
start, the camera was started and the stopper 
pulled. This released the oil at a fairly fast 
rate. In order to determine the amount released, the 
fluid level in the cannister was measured before 
and after the release. 


The underwater experiments were considerably 
more difficult. The experiment had to be set up in 
such a way that no air pockets remained in the tank 
and no kerosene escaped through the feed lines 
connecting the tank to the cannister. To do this, 
water was first placed in the cannister to a level 
that would prevent the kerosene from getting into 
the feed lines. Then, enough kerosene, which was 
dyed using enamel paint, is added to fill the 
cannister. The rough surface was then placed on 
the tank and the tank was almost completely filled 
and then tilted to allow the air to escape from the 
pockets. The filling of the tank was completed in 
this position. When all the air escaped out the 
top and the tank was full, it was rotated back to a 


horizontal position and leveled. Since there was 
sufficient pressure in the tank to cause the surface 
to deflect at the center, water was bled off in an 
attempt to ensure that the surface was flat. Some 
degree of error was probably caused by this deflec- 
tion. Once the tank was leveled and the pressures 
equalized, the experiments proceeded in the same 
way as the over ice experiments. 


To check for the presence of viscous effects, 
experiments were run with glycerin-water and 
kerosene-oil mixtures. These experiments are 
described in the Appendix. 


Lil. Lneor 


The theory that governs the spreading was 
developed as a result of data from the experiments 
described in the last section and the available 
field data. 


In determining what will happen to oil spilled 
in the Arctic, two cases must be considered. In 
the case of a pipeline breakage, the oil will spread 
over ice and tundra. In the case of a supertanker 
accident, the oil may spread under the ice. This 
happens because Arctic ice has a specific gravity 
of about .8 and North Slope crude has a specific 
gravity of about .9 . Therefore, at an ice-water 
interface the oil will, on hydrostatic considera- 
tions, tend to spread under the ice. The spreading 
in both cases follows the same theory and the re- 
sults differ only by a factor determined by the den- 
sity ratio of water and oil. 


Patureau's theory states that pockets_in the 
ice are filled to the roughness height, <h > , 
and it was shown experimentally that the value of 
the viscosity is not important - This means that 
as pockets in the ice are filled additional fluid 
flows over the filled pockets frictionlessly and 
the drag can act only at the edge over an area pro- 
portional to the radius multiplied by the roughness 
height. When attempting to correlate the labora- 
tory data for over and under ice and the field data, 
it was found that the data would nepoycale with the 
volume released. It was found that the rate of 
spread for a given surface could instead be corre- 
lated with using the rate of oil release, Q . 
Since the size of the final pool depends upon vari- 
ous statistical parameters, one can see that for a 
given volume and surface there will be an average 
pool size and an associated variance. In order £9) 
determine the accuracy of the averages, Patureau 
developed an expression relating the standard devia- 
tion of the size with the average size. 


Fig. 3 shows a sketch of oil being released 
on a rough surface with a roughness height <h > . 
The driving force is gravity, which causes a pres- 


sure at the outside edge of the spill. Thus, 
Be ~ pA ~ (pAgh)hr (1) 
Since 
en (2) 
r 


the gravitational force from Eq. (1) can be 
expressed as 


2 
F ~ pig 0 : (3) 


r 


The pressure drop caused at the outer edge as the 
oil flows will be 


2 
d, ~ Ce ov (4) 


where C is a friction coefficient of order one. 
The area over which the drag acts if the roughness 
height, <h >, times r . The retarding force be- 
comes 


F~ our <h > (5) 


Replacing U in Eq. (5) with r/t and setting 
Eq. (5) equal to Eq. (3) gives 


2 3 
pdg ED. DE n> (6) 
18 t 
This result can be written as 
1/6 
2 
ety BEOE 22/3 (7) 
<h > 


Eq. (7) can be non-dimensionalized by defining a 
scaling length 


* = 
Taee=— <he> (8a) 
and a scaling time 
1/4 
* <h >! 
t = (sz) (8b) 
2Q 


To obtain an expression for the final size and 
the associated ee toe deviation several assump- 
tions must be made The first of these assump- 
tions equates the average of a sample from the sur- 
face with the ensemble average over the whole sur- 
face. It is assumed that the elevation of points 
on the surface are normally distributed about some 
mean. This assumption is in only fair agreement 
with the experimental surface, which actually is 
closer to a lag-normal distribution fapd good 
agreement with the real ice surface) _Under 
this assumption, the roughness height, <h > , and 


the variance of the height, oy » can be related 


by the simple equation: 


il // 22 
<h > = (2 ) (9) 
27 


Thirdly, it is assumed that pockets in the ice 
surface are filled to the mean roughness height (see 
Fig. 4). This greatly simplifies the theory, but 
is only a crude approximation. The accuracy of this 
assumption can only be determined by experiment. 
The fourth main assumption is that the number of 
pockets filled is always large enough to give a 
statistically meaningful result. This assumption 
is realistic. For yisblicity the pockets have a 
rectangular shape (square in the case of an 
isotropic surface). This was done strictly for 
computational reasons. For the purpose of this 


paper the surface will be assumed isotropic and the 
spill will be circular. The isotropy assumption 
will be shown to be a good assumption by comparing 
the power spectra along perpendicular paths of the 
experimental surface. Lastly, it was assumed that 
a two-dimensional power spectrum could be formed 

as the product of two-dimensional spectra. This 
assumption has no mathematical justification, but 
even if it were poor, no real error in the results 
would appear. 


The final results of the theory are expressed 
in terms of parameters which are derived from the 
statistics of the sea ice surface. The first of 
these is the variance of the roughness, 0, . Next 
are two parameters k and T derived from the 
power spectrum of the surface. For a one-dimension- 
al spectrum, T is the wavelength of the peak of 


the spectrum and k is determined from the follow- 
ing equation: 
2710. 
Lo h -2T 27 
PC) Sebati /D: exp ( i/? cosh(~775) (10) 
k Tk Tk 
To simulate sea ice, the product kr? should be 


greater than 50, or, in the two-dimensional case 


(kT?) should be greater than 2500. From (kt?) ; 
an additional parameter, 9, can be determined 
from Fig. 6. 


Following the assumption that the oil fills 
the pockets to the roughness height on the average, 
we get the equation 


<Ysea)<S e> <he>. (11) 
max 
or, for a normal distribution, 
oO, 1/2 
< = > 7 
Ve Shae 4, Se 


from Eq. (8). 


To non-dimensionalize the results, Patureau 
introduced the,concept of the ‘most probable 
pete) 
pocket' size, » that is, the most common size 
for a pocket on the ice. Defining the surface area 
of a most probable pocket as 


2 
A THeay (13) 
and the volume of a most probable pocket as 
2 (Ae) yi/2 
a, eye a 
abe mkar 26 ary i’) 
P T 


we find the following dimensionless equation from 


Eiqem@z)) hc 
max” _ p> eee 
Sip 73/2 vas 


To provide a measure of the accuracy of the 
above equation, an expression for the variance of 
the surface area, a » was derived. This 


(15) 


expression is obtained by considering the pockets 
in the ice to be filled sequentially and independ- 
ently. For ease of understanding, the derivation 
will be for the one-dimensional case. This can 
then be easily extended to the two-dimensional case. 


Since the pockets each have a volume variance, 


fo} 


Vv » and a length variance Ps and since they 


mp 
are filled independently, the total variance oy 


will be expressed as the sum of the variance of the 
pockets filled: 


(16) 


The total length of the average one-dimensional 


Spice OK > is equal to the length of a mean 
pocket, X » multiplied by the number of pockets 
filled,N. 
<x NX (17) 
max mp 
oy (5 o, (18) 
mp 
Using Eq. (16) 
O, = 0k No, (19) 
mp mp 
or 
1/2 =) 1/2 ey? (19a) 
10) = N oO 
x Vv Vv 
mp mp 


In dimensionless form, this becomes, using Eq. (16): 


1/2 1/2 
eu <x S °y 
x = max é mp (20) 
xX V 
mp mp mp 
In the two-dimensional case, one finds 
1/2 
0 1/2 Bot > 
s x 
ge ce Oils rer (21) 
mp mp 


Where the proportionality constant can be defined 


as 
oy 1/2 
(ee ep 
mp 


with © as shown in Fig. 5. 


(22) 


It is seen that 9 


depends only upon the product kT . This yields 
Efe) 
Oo 1/2 ase > 
s ~ x 
oe = 0 5 (23) 
mp 


which is the result presented by Patureau. 


IV. Results 


Figs. 6 and 7 show the power spectra for the 
surfaces 1 and 2, used in the experiments. The 
isotropy assumption can be checked by comparing the 
spectral points determined from paths taken per- 
pendicularly. The two directions are distinguished 
by the use of open and closed symbols. By noting 
that there is little difference between the two 
directions, one can see that the isotropy assumption 
is a good one. Also shown on these figures are 
on eee Landen Om. 


Figs. 8 and 9 show that the data for the final 
pool size for surfaces 1 and 2 respectively is in 
good agreement with Patureau's theory. The values 
of the scaling area and scaling volume are shown on 
the figures, where S and Vy) are as calculated 
in Eqs. (13) and (14)™P The peta lines on the 
graphs are the average values of the final size 
from Eq. (15). The dashed lines on either side 
represent one standard deviation on each side of 
the average and were determined from Eq. (23)... 3.00 
these graphs, the open symbols represent over ice 
runs and the filled symbols indicate under ice runs. 
It appears that the under ice pools are slightly 
smaller. This effect is within experimental error. 


Fig. 10 shows the dimensionless radius versus 
dimensionless time for both over and under ice 
cases using the results of Section III. Also 
shown gn this graph ts field data from Glaeser and 


Vance and McMinn The equation that best 
describes the line in Fig. 11 is: 
1/6 
Se yey) 


The under ice data shown in Fig. 13 is plotted 
after allowing for a starting time. It was found 
that in the under ice case there was a significant 
lag between the time the stopper was pulled and the 
actual beginning of the spread. This was due to 
the time for the kerosene to float up through the 
water to the surface. This starting time was found 
to be about 1.5 seconds. 


For the field data, the roughness heights are 
not known. Thus, the field data was fit to the 
solid line in Fig. 13 by adjusting the value of the 
roughness height. In most cases, the value of 
roughness height appears consistent with visual 
observations made in the field. 


To expand the results of the laboratory 
experiments, one must know the characteristics, 
particularly the roughness height, of the Arctic 
ice. From laser profilometer traces taken in areas 
of high spill potential it has been determined that 
the roughness height varies greatly but is seldom 
less than three centimeters or greater than 60 
centimeters The size of a spill would be on the 
order of 10° cubic meters. Fig. 14 shows how the 
final radius of a spill varies with the roughness 
height and the total volume released. For an 
average volume and roughness height the final radius 
will be on the order of 300 meters. 


Fig. 


Fig. 


Fig. 


Fig. 
Fig. 
Fig. 
Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


BIBLIOGRAPHY 


Blackman, Ralph B., and Tukey, J. W., The 
Measurement of Power Spectra, Dover, N.Y.,1959. 


CRREL, SEC Arctic Environment Data Package, 
Corps of Engineers, Hanover, N.H., 1970. 


Glaeser, John L., and Vance, George P., A Study 


of the Behavior of Oil Spills in the Arctic, 
U.S. Coast Guard, Washington, D.C., 1971. 


Granger, C.W.J., Spectral Analysis of Economic 
Time Series, Princeton University Press, 


Princeton, N.J., 1964. 


Hoult, David P., ed., Oil on the Sea, Plenum 
Press, N.Y., 1969. 


McMinn, T.J., Crude Oil Behavior on Arctic 
Winter Ice, U.S. Coast Guard, Washington, D.C., 
1972). 


Patureau, Jean-Pierre, Statistical Approach 


for Determining the Extent of an Oil Spill Over 
a Rough Surface, M.I.T. Master's Thesis, 1972. 


Suchon, Walter, An Experimental Investigation 
of Oil Spreading Over Water, M.I.T. Master's 


Thesis, 1970. 


O'Dea, Stephen R., Experimental Determination 


of the Extent of an Oil Spill on Arctic Ice, 
M.I.T. Master's Thesis, 1973, Dept. of 


Mechanical Engineering. 
LIST OF FIGURES 


1k Schematic of Spreading Apparatus (over 
Ice Spread) 


2 Schematic of Apparatus to Measure Surface 


3 Section of Oil Spreading Across a Rough 


Surface 
4 Average Longitudinal Dimension of a Pocket 
5 Q vs Dimensionless Group ioe 
6 Power Spectrum vs.Wavelength for Surface 1 
7 Power Spectrum vs.Wavelength for Surface 2 
8 Dimensionless Volume Released vs.Surface 
Area for Surface 
9 Dimensionless Volume Released vs. Surface 


Area for Surface 2 


10 Dimensionless Radius vs. Time for Field 
and Laboratory Data 


ll Average Areal Coverage vs. Volume 
Released for a Spill on Arctic Ice 


<h > 


LIST OF SYMBOLS 


roughness height of a surface 


parameter describing rough surface that 
is determined from power spectrum 


number of pockets filled 
rate of oil release 


average value of the surface area of a 
spill 


surface area of a mean pocket 
wavelength of peak power spectrum 
total volume of oil released 
volume of a mean pocket 


average value of the length of a spill 
(one-dimensional case) 


length of a mean pocket (one-dimensional 
case) 


roughness height of a surface 


density ratio of water and oil, 





density of oil 

density of water 

variance of the roughness of a surface 
variance of the surface area for a spill 


variance of the volume for a spill 
variance of the volume of a mean pocket 


variance of the length of a spill (one- 
dimensional case) 


non-dimensional parameter which is the 


1/2 


slope of the o VS. <S > relation 
s max 







MIRROR 


CANNISTER 


CAMERA 





ICE 





SIMULATED 


FIG. 1 SCHEMATIC OF SPREADING APPARATUS (OVER 


ICE SPREAD ) 


FILLED 
POCKET 


RELEASE 
RATE =Q 






DRAG ONLY AT 
LEADING EDGE 


FIG. 2 SECTIONOF OIL SPREADING ACROSS A 
ROUGH SURFACE 


(ONLY 8 OUT OF 2a SHOWN FOR EASE ) 


oe SPOCKE? 


hex) 





¢h> =0 





EXPECTED VALUE a, 


FIG. 4 AVERAGE LONGITUDINAL DIMENSION OF 


A POCKET 


. 30 


lls) 





File, 2 





SURFACE 


SCHEMATIC OF APPARATUS TO MEASURE 
SURFACE 
104 10° 10° . 
(KT? 


FIG.5 @ vs DIMENSIONLESS GRAND 1x2 2 


if 


10 








A (in.) 


FIG.6 POWER SPECTRUM vs WAVELENGTH FOR 











SURFACE 1 
2 
10 = 
SYMBOL RUN ' 
i) al 
ra 2 
° 3 
o 4 
A 5 
@ 6 
2 Foe THEORY 
10 Wie EQ= (15) pal 
a 
Sop 











Af 
V 
mp 


FIG. 8 DIMENSIONLESS VOLUME RELEASED vs SURFACE 


AREA FOR SURFACE 1 





16 8 4 Z 1 
7 (tits) 


FIG. ‘7 POWER SPECTRUM vs WAVELENGTH FOR 


SURFACE 2 























10° [ | 
S 
Cc SYMBOL RUN 
mp a a 
A 8 
@ 9 
0 15. 
a 16 
O° iy) 
iQ) = We 7 
eee 
THEORY 
EQ. (15) 
L | 
1 10 10 
ale 
V 
mp 
FIG. 7 DIMENSIONLESS VOLUME RELEASED vs SURFACE 


AREA FOR SURFACE 2 


<h> 





i0 
SYMBOL 
0 
@ 
lo” ‘ 
Fadi 
=f 
fi 
| 10 


FIG. 10 DIMENSIONLESS RADIUS vs TIME 


(meters) 





| ud ot 
8 
ao 
RUN q 
e 
1-4 
5-6 
14-16 
7-9, 17-18 
8 
+a @) eo 
$A 
sen 
B, ae 
e 
<A SYMBOL <h >(m) Q(m°/sec) 
a 
KS A .010 .0049 
4 038 0049 
a 045 0049 
A 085 0049 
a . 00030 .057 
N 00030 .113 
a .00036 062 
g 0042 .0079 
a 025 .00048 
lo” 10? Io" lo” 


4 

5/7 VV. 
(<tyK4q2)4 
FOR LABORATORY AND FIELD DATA 


Z = EFFECTIVE 
ROUGHNESS HEIGHT (m) 


FIG. I} AVERAGE AREAL COVERAGE vs VOLUME RELEASED 


FOR A SPILL ON ARCTIC ICE 


-- NOTES -- 


DATE 
5 DUE 





Pam: 502.55:665.6 | HOU 
ee UU 









_ HOULT, David P. 


The spread of oil in the Arctic 







#13669 


Boreal Institute for Northern 
Studies Library 

CW 401 Bio Sci Bldg 

The University of Alberta 

Edmonton, AB Canada T6G 2E9 








niversity of Alberta Library 


wii 


0 1620 0329 2081 


